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Abstract: The purpose of this talk is to sketch some recent progress which has been made in calcu- 
lating non-perturbatively the reflection factors for the sinh-Gordon model restricted to a half-line by 
integrable boundary conditions. The essential idea is to calculate the energy spectrum of boundary 
breathers in two independent ways; firstly by using the boundary bootstrap and secondly by quan- 
tizing the classical solutions corresponding to boundary breathers. Comparing these two calculations 
provides a way to determine the dependence of the reflection factors on the parameters introduced 
into the Lagrangian by the boundary conditions. The basic idea is illustrated using a massive free 
scalar field with a linear boundary condition confining it to a half-line. 
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1. Introduction 

Although boundary conditions are extremely im- 
portant for field theory — think of monopoles or 
vortices, for example, or instantons — the study 
of integrable quantum field theory restricted to 
a half-line or an interval by integrable boundary 
conditions is not yet fully developed. In particu- 
lar, many interesting areas are yet to be explored 
and the foremost of these is the classification of 
the reflection factors which determine the man- 
ner in which a particle rebounds from a bound- 
ary. In 1993 Ghoshal and Zamolodchikov @ 
ignited activity in this subject by analysing the 
possible behaviour of the sine-Gordon model re- 
stricted to a half-line. They discovered a two- 
parameter family of boundary conditions which 
preserved integrability, in the sense of preserv- 
ing after suitable modification all the energy-like 
conserved quantities, while violating those which 
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were momentum-like. Clearly, the same facts 
would be true of the sinh-Gordon model, which 
is even simpler in the sense of having a particle 
spectrum consisting of a single scalar particle and 
nothing else. The reflection factors for the states 
belonging to the sine-Gordon particle spectrum 
(solitons, anti-solitons and breathers) were com- 
puted although the relationship between two pa- 
rameters appearing in the reflection factors and 
the two parameters appearing in the classical La- 
grangian specifying the boundary conditions re- 
mained largely obscure. The purpose of develop- 
ing the techniques reviewed here was to remove 
this obscurity. It seems such a simple matter that 
it is surprising that it has taken so many years 
to find a way to resolve it. 

Apart from the troubles relating Lagrangian 
parameters to the parameters appearing in reflec- 
tion factors there are other mysteries too. For ex- 
ample, the sinh-Gordon model is the simplest of 
the affine Toda field theories and yet it is the only 
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one in the ade series of cases with a set of bound- 
ary conditions with free parameters. For all the 
rest, the possible integrability-preserving bound- 
ary conditions unexpectedly form a discrete set 
H |l0[ ^ . It is not at all clear what the reflection 
factors are which correspond to each member of 
this discrete set, although a variety of reflection 
factors have been found using algebraic means 
(see for example [0). Our long-term goal is 
to find a complete classification and a full un- 
derstanding of the curious pattern of integrable 
boundary conditions discovered in . 

For a review of the background to some of 
these questions see 

2. A free scalar field with a linear 
boundary condition 

In this section, to set the scene, consider a free 
massive scalar field confined to the left half-line 
(x < 0) by a linear boundary condition at x = 0. 
Its equation of motion and boundary condition 
are: 

(d 2 + m 2 )(j) = 0, x < 0; d x (j) = -\<p, x = 0. 

(2.1) 



Clearly, there is a set of solutions to (2.1) cor 



responding to a superposition of left and right- 
moving waves of the form 



(e lkx +Ro(k)e- lkx ) +c.c. (2.2) 
in" + k 2 and a 'reflection factor' given 



by 



Rq 



ik + A 



(2.3) 



ik — A 

If A < there is also a periodic solution of the 
form (whose existence is indicated by a pole in 
Ro at k = —iX), 



= Acosuite Xx , u> z — m 2 - A 2 , 



(2.4) 



provided also > A > — m. Clearly, this solution 
is slightly delicate and fails to exist at all for a 
massless free field. It is a harmonic oscillator 
'glued' to the boundary. Its energy is computed 
to be £ where 



£ = 



dx (j) 2 + <j)' 2 + m 2 2 ) + A0 2 (O, t) 



lo 2 A 2 
2A 



Quantizing this periodic system leads to a tower 
of states with energies (n + l/2)u>. 

On the other hand, in the context of inte- 
grable field theory, we know that the S-matrix 
for the free massive particle is unity and the re- 
flection factor is identical to the classical ver- 
sion ( |2.3| ). That is, setting u> — mcosh.0, k = 
msinh# we have the relationship between 'in' 
and 'out' states 



\B >out=Ra(k)\ -9 >i 



(2.5) 



where 'in' refers to a particle of rapidity ap- 
proaching the boundary and 'out' refers to the 
same particle at a later time after reflecting from 
the boundary. The pole in Rq at k = —i\ indi- 
cates a boundary bound state at an imaginary ra- 
pidity (0 — iijj — — i sin^ 1 (A/m)) and the bound- 
ary bootstrap [|l8| [ij], ^, indicates that the en- 
ergy of the excited boundary is given by 



Si = £q + m cos ip = £q 



(2.6) 



with an associated reflection factor again equal 
to i?o- The reason for this is that the reflection 
factor i?i for the particle rebounding from the 
excited boundary is given by 

Ri{6) = S(0 - itlj)R a {6)S(9 + itp) (2.7) 

and the S-matrix is unity. Repeating the process 
leads to the tower of harmonic oscillator states. 

So, the basic idea which we wished to exploit 
was to follow this line of reasoning in the sinh- 
Gordon model, constructing a tower of bound- 
ary states in two different ways: firstly, using the 
formal bootstrap, and secondly, using a direct 
quantization of periodic solutions corresponding 
to an 'oscillator' attached to the boundary. It is 
expected that the same approach will be appli- 
cable to the other Toda theories as well. 

3. The sinh-Gordon model 

We begin by establishing notation and reminding 
ourselves of certain relevant facts about the sinh- 
Gordon model on a half-line. Its field equation 
and integrable boundary conditions are 

d 2 4> = -^S^ sinh(V2/30), x < (3.1) 



P 2 



eie 



)■ 



0. 
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where m is a mass parameter (which we shall set 
to unity), (3 is the bulk coupling, and eo and e\ 
are the two parameters introduced by the bound- 
ary conditions. Notice that when eo ^ ei the 
bulk symmetry under the transformation — > 
— is lost. Although it is not strictly necessary 
to assume that | e.^ | < 1, in what follows it will 
often be useful to put 

€i = cosaiir, i = 0, 1, < ai < ao < 1. (3.2) 

The basic question is to give a complete charac- 
terization of the reflection factors as functions of 
the bulk coupling and the two boundary param- 
eters. 

The sinh-Gordon model is not free but its 
S-matrix is well-known and may be written com- 
pactly as follows 

1 „,„s P 2 /2ir 



4. Using the boundary bootstrap 



S(6 12 ) = - 



(B)(2-B 



-, B(J3) 



1 



-/? 2 /47T' 

(3.3) 



(3.4) 



using the block notation, 

_ sinh (0 12 /2 + i'Kx/4) 
^ ~ sinh (6 12 /2 - iirx/A) ' 

introduced in ||]. The S-matrix is invariant un- 
der the transformation (3 — > &n/f3. 

The reflection factors may be deduced from 
those computed by Ghoshal for the sine- 
Gordon model, and have the form 

(l)(l + S/2)(2-fl/2) 



R(0) = 



(3.5) 



{l-E)(l + E)(l-F)(l + F)' 

where E, F are two functions which are inde- 
pendent of the rapidity but depend upon the 
bulk coupling and the boundary coupling param- 
eters in a manner which needs to be determined. 
The classical limit of (3.5) as j3 — > has been 
calculated independently [jl0| and implies E — > 
ao + ax, F — ► ao — Oi, using the notation (|3.2|). 



Indeed, (|3.5|) is the simplest guess having the ap- 
propriate classical limit and satisfying the 'cross- 
ing unitarity relation' [[18], |l6[ 

R(0 + in/2)R(6-iir/2)S(2O) = 1. (3.6) 

We are also intrigued to know what happens to 
the weak-strong coupling invariance enjoyed by 
the S-matrix factor. 



The expression (3.5) has also been checked 
perturbatively in various ways J|(| [?], [|] but 
the details will not be described here. 



Examining the expression (3.5) reveals that it 
may have poles at 9 = itp, < ip < ir/2 and 
that these poles cannot come from the numer- 
ator factors, only from zeroes of the denomina- 
tor. Moreover, the existence of suitable poles de- 
pends on the value of E or F. For example, with 
1 < E < 2 there are poles at = i(E - l)ir/2. 
What is more, knowing the classical limit of E 
and F is useful because it informs us that it is 
plausible, given the ranges of a and Oi, that 
E and F do not lie simultaneously in the range 
(1, 2). This fact is true for the classical limit, and 
the order 1 perturbative calculations do not up- 
set it. So, we shall proceed as though it were true 
and see that the resulting expressions for E and 
F are consistent with it. For the symmetric case, 
(ao = ai), F = 0, and the question does not even 
arise. More details of these calculations may be 
found in g and Jl§. 

Then, provided 1 < E < 2 we may compute 
a tower of states using the boundary bootstrap 
(2.7). The result is a sequence of states whose 



energies are given by 



S n+l = S n + m cos ~(nB- E + l), (4.1) 



and for which the associated reflection factors 
have the form 



Rn (0) — Rd 



{l + E + B){l-E~B) 
(1- E + nB){\ + E - nB) 
1 

(1- E + (n~ 1)B)(1 + E-(n- 1)B) ' 

(4.2) 



where 



Rl 



(!)(! + B/2)(2- B/2) 
(1-F)(1 + F) 



(4.3) 



is n-independent and never contains contributing 
poles. 

Notice that the number of new types of pole 
is limited because, even if E is initially within 
the appropriate range, eventually E — 1 — nB 
will move out of range. It seems natural to in- 
terpret the pole at i(E — 1 — (n — l)B)ir/2 as a 
'crossed' process in which the state with energy 
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£ n drops down to the state at energy £ n _i. If 
we interpreted this pole as a new state then the 
reflection factor corresponding to it would con- 
tain higher order poles and these in turn would 
need to be explained. 

The next part of the argument requires a set 
of classical solutions to the system of equations 



(3.1) 



5. Classical boundary breathers 

The classical periodic solutions which we shall 
use may all be described within Hirota's ansatz. 
That is, for our purposes it is convenient to set 



\/2. r 

~T ln v 



(5.1) 



where 



T a = 1 + {-) a {Ei +E 2 + E 3 ) + A 12 E X E 2 

+Ai 3 EiE 3 + ^23^2-^3 + Ai 2 Ai 3 A 2 3EiE 2 E 3 , 

(5.2) 



and 



iTi a p x-\-b p t-\-c p 

n/p — e 

a p = 2 cosh p p , b p — 2 sinh p, p 



(5.3) 



At this stage, the constants c p are not deter- 
mined. 

On the whole line, these solutions are not 
particularly interesting because for each t they 
develop a singularity at some x. Moreover, ow- 
ing to the singularities their total energy will be 
infinite. However, we are interested in periodic 
solutions defined only on the region x < and 
therefore it is possible that solutions satisfying 



1 However, there is a bit of a mystery, or at least an 
inconsistency of interpretation between this case and the 
free field case. Here, we have in mind that there ought 
to be a finite number of boundary states just as there 
are a finite number of sine-Gordon breathers at most, at 
a given bulk coupling (see J22| for a review); only in the 
classical limit do we expect to have infinitely many bound 
states. On the other hand, it seems natural to have the 
infinite tower when we are discussing the free field bound- 
ary states. Thanks to Patrick Dorey for a discussion on 
this issue. 



the boundary conditions can be non-singular ex- 
cept at points beyond the boundary; in effect, the 
singularities are hidden away in the region x > 0. 
This hope turns out to be realised rather neatly 
and the details are given more fully in The 
trick is to select E\ and E 2 to be time depen- 
dent and periodic (ie, inevitably complex), and 
to satisfy E\ — E 2 , while taking E 3 to be real. 



Then, the solution given via (54) is real and the 
boundary conditions can be satisfied. 

For the particular case where eo — e\ = e, it 
is enough to set E 3 = and we have, 



T a = 1 



(_). 



} 2x cos p 



tanp 

e 4x cos p e + cos P 
e — cos p 



cos(2i sin p) 



cos p 



cos p 
(5.4) 



The expression (5.4) is clearly periodic, with pe- 



riod 7r/ sin p. However, for the corresponding so- 



lution constructed via (5.1) to be non-singular, e 
and p must be constrained as follows: 



T < e < 0, cos p < —e. 



(5.5) 



As expected, e should be negative (meaning that 
there is a competition between the different signs 
of the boundary and bulk contributions to the 
energy), but not too negative (cf the free field 
case). More of a surprise is the constraint on p 
since it means that there is a minimum frequency 
for the breather but, at that minimum frequency 
the amplitude of the breather has shrunk to zero. 
This behaviour is different to the behaviour of 
the sine-Gordon breathers since for those the fre- 
quency drops to zero as the amplitude collapses 
to zero. It seems the boundary breathers are 
still breathing faintly despite their collapse when 
cos p = — e. This behaviour is reminiscent of 
a standard oscillator whose amplitude may be 
tuned to zero independently of its frequency. In 
a nonlinear system for a given boundary condi- 
tion the frequency governs the amplitude of the 
periodic solution. 

In the general case, the solutions are sub- 
stantially more complicated. The boundary con- 
ditions require 



2 " 



S q- 
tan p 



1 1 - q 2 q\ 1 — q 2 q 2 _ 



j 2 1 — q~ 2 q\ 1 — q 2 q 2 _ 



(5.6) 
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where 

p 17 

q = tan-, q± = tan-(a ± a x ). (5.7) 

The general boundary breathers have no singu- 
larities in the region x < provided 

7T 7T 

cos — (clq + a± ) < 0, cos — (a — ai) > 

7T 

< cos p < — cos — (do + tii). (5.8) 



The restrictions f p-aj) were found numerically us- 
ing Maple assuming that the parameters ao and 



ai lie in the range indicated in (3.2). The energy 
of the general periodic solution given by the data 



( 5J3 ) is given by 
£ b = 



1 na .. . , 

^ ^1 + cosp- -(sm— +sm— ) 



770,1 ^2 

(5.9) 



6. Semi-classical quantization 

Now that we have a collection of periodic solu- 
tions the next step is to quantize them to deter- 
mine their energy spectrum. We shall do this by 
adapting the WKB methods proposed long ago 
by Dashen, Hasslacher and Neveu p3| . Here it 
will be enough to outline the prescription and 
more of the details are to be found in || [l2|] . 
The procedure is certainly non-perturbative but 
we have no way of telling at present if it is also 
exact. In the case of the sine-Gordon model in 
the bulk, the WKB methods do give exact results 
in agreement with calculations made within the 
eight vertex model (for example, see 21 ). It re- 



mains to be seen whether this feature will persist 
in the presence of boundaries. As we shall see, 
there are several surprises which render the anal- 
ysis for the general periodic solution manageable. 

The first ingredient we need is the classi- 
cal action integrate over a single period of the 
boundary breathers. Since we know the energy 
already, it will be necessary to calculate the ap- 
propriate integral of 4> 2 . Thus, 

r T ,0 

Sdass = / dt I dx(t> 2 - TZ h 



8-7T / a + ai 



— TEh 



(6.1) 



This is a remarkable result because the first part 
of the expression depends only on the combina- 
tion clo + cli. For the symmetric case ao = oi, the 
expression for the action is relatively straightfor- 
ward to compute analytically by a suitable se- 
quence of changes of integration variables. How- 
ever, in the general case the result was deduced 
numerically and the rather convincing numerical 
evidence we have available is summarised in the 
appendix to jl2| . It would be nice to have an an- 
alytic derivation of fl6.1|), however it has eluded 
us so far. 

Besides the classical action we shall also need 
the linear perturbations around the breathers re- 
garded as a background. In other words put 
4> = <po + 7], where tfio is a breather and solve 
the following linear equations for j]\ 

d 2 r] - d 2 x ri + At] cosh(V2/3(/) ) = 0, x < 



d x rj + (e Q e~ 



■04>o/V2 . 



eie Wo/v^j 



0, x = 0. 
(6.2) 



The solutions in which we are interested are asymp- 
totically plane- wave solutions, which means that 
as x — > — oo they have the form, 

77 - e~ wt (e lkx + R{k)e' lkx ) , u 2 = k 2 + 4, 

(6.3) 

where R is the classical reflection factor in the 
breather background. Actually, we shall need the 
classical reflection factor in the classical ground 
state background as well (this particular reflec- 
tion factor, incidentally, is the classical limit of 



Ghoshal's formula (3.5). Fortunately, a further 
use of Hirota's method, adding two infinitesimal 



exponentials to (5.2) in the standard manner, al- 
lows a straightforward derivation of the reflection 
factors we need. They are, 



Rb-c 



cathcr 



ik-2 



COS/9 



ik 4 
ik 



ik-2 



2 cosp J ik 
- 2 cos a + ik 



2 

2 cosa_ 



R 



ik 



ik + 2 cos a+ ik — 2 cos a_ 
-2 ik + 2 cos a + ik + 2 cos a_ 



ground — 



ik + 2 ik — 2 cos a + ik — 2 cos a_ 



(6.4) 



and, for notational convenience it is useful to de- 
fine a± = 7r(ao ± a\)/2. The second of the re- 
flection factors given in ( |6.4|) allows us to deduce 
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the classical limit for E or F quoted earlier. An- 
other surprising feature which we did not expect 
is that the ratio of the two classical reflection fac- 
tors does not depend on the combination ao — a\ . 

The linear perturbations of the breathers do 
not have the same periodicity as the breathers 
themselves. Rather, their frequencies provide the 
modification to the action which is to be used 
in the WKB aproximation. Thus, we need to 
calculate a quantity A defined by 



A ^ y breather Aground) ; 



(6.5) 



and, in terms of this the quantum action is de- 
fined to be 



Sq u — S<L 



lass 



A. 



(6.6) 



The evaluation of A is not quite straightforward 
although we are fortunate to have the old work of 
Dashen, Hasslacher and Neveu as a guide. The 
method we adopted runs as follows. First dis- 
cretize the sum over the frequencies and then 
regulate it by removing an infinite piece. It is 
convenient to use a Dirichlet boundary condition 
at x = —L, so that t](—L, t) — 0, effectively plac- 
ing the system in a box without changing the 
boundary condition at x = 0. Subsequently, the 
limit L — ► oo will be taken at the end of the calcu- 
lation. Removing the infinite part of the discrete 
sum corresponds to making a choice of normal- 
ordering in a perturbative approach to the field 
theory. 

Once the Dirichlet condition is imposed the 
possible discrete values of the parameters k ap- 
pearing in ( |3.3| ) are given by 



-2iLk £ 



-2iLko 



flo(*o), (6.7) 



where the subscripts B and refer to the breather 



and ground state, respectively. Typically, (6.7) 
have infinitely many real solutions together with 
a small number of purely imaginary solutions. 
The latter do play a role which will not be dis- 
cussed here. The details may be found in fl2]| . 
For large L, it is useful to note that the differ- 
ence between ks and fco is small, prompting us 



to write for the solutions to (6.7) 



(fcn)n - (k )n + T «((fco)n), (6.8) 



where the function k satisfies 

ik + 2 cos p ik + 2 cos o_| 



2i K (fc) _ 



ik - 



(6.9) 



2 cos p ik — 2 cos a_| 

Again, we notice the curious and surprising fact 
that only the combination ao + a± appears. 

The detailed calculation of A may be found 
in H p^|; here it is enough to quote the result 
which is pleasingly simple: 

2 



(cos p - 



p sin p 



+ {a + - 7r/2)sina + ) . (6.10) 

The WKB prescription instructs us to define 
£ qu = —dS qu /dT, and then set 



Squ + T£„ u = 2nir 



(6.11) 



where n is a (positive) integer, or zero. In other 
words, assembling all the ingredients, we have 
the quantization rule 

2mr=i(p-7r/2) + ^( a+ -7r/2), (6.12) 

from which we deduce a set of angles p n , and 
thence a tower of energies £ n satisfying 



Cn+i = tn H „ sin — r x 

cos - y-p{a + - 7T/2) - (n 



1/2)5 



(6.13) 



which should be compared with (4.1). Before 
making the comparison it is perhaps worth not- 
ing that as (3 — > 



Pn 



(tt - a+), £ n -^(n + l/2)u , (6.14) 



where ujq = 2 sin a + is the lowest possible breather 
frequency. 

Comparing (|6.13| ) with (4.1) for several val- 
ues of n yields two pieces of information. First 
of all the mass of the sinh-Gordon particle must 
be given by 



(6.15) 



and, secondly, the parameter E appearing in ( |3.5| ) 
for the reflection factor must be related to the 
boundary parameters and the bulk coupling by 

E{fi, a , oi) = (a + oi)(l - B/2). (6.16) 
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This expression agrees perfectly with the classi- 
cal limit as (3 — > 0, and with the result given 
by Ghoshal and Zamolodchikov for the Neumann 
boundary condition. There, ao = di = 1/2 and 
E = 1 — B/2 as they claimed. The expression 
fl6T6| ) also agrees with one loop perturbative cal- 
culations , and with the known results at 
the special point B = — 2 Q (see @). 

As far as the other parameter F is concerned 
we know its classical limit, and we also know its 
expansion to order 1 , at least to first order in the 
difference ao — a\ (see 0). Unfortunately, when 
the bulk Z2 symmetry is broken the perturba- 
tion calculations are substantially more compli- 
cated, and indeed unfinished for arbitrary ao,ai 
(see ||). These facts are, however, consistent 
with the expression 

F(J3, a , oi) = (a - oi)(l - B/2). (6.17) 

In any case, had we chosen to replace 01, say, by 
— ax, we would expect E and F to interchange 
(since E by itself is not invariant under such a 
change) . 

Notice that since < 1 — B/2 < 1, the as- 
sumptions we made concerning the regions in the 
parameter space for which the reflection factor 



If we define a new triple by making the change 



47T . 
/J 3 



(fl, ao, ax) -> ^(J3, a , ax), (7.1) 



(3.5) has poles are vindicated. 

One slight criticism of the technique is the 
requirement on the boundary parameters to lie 
in the region where the boundary breathers ex- 
ist. Otherwise, strictly-speaking we could de- 
duce nothing about the reflection factor since 
the comparison we have utilised would not be 
valid. However, it is perhaps worth noting that 
in the limit where ao = ax — > 1/2, the result does 
agree with the earlier conjectures of Ghoshal and 
Zamolodchikov concerning the Neumann bound- 
ary conditions despite the obvious fact that the 
Neumann boundary condition does not support 
boundary bound states. 

7. Weak-strong coupling duality 



We mentioned that the S-matrix describing the 
scattering of the sinh-Gordon particle is invariant 
under the exchange f3 — > 4tt//3. It is interesting 
to notice that there is a sense in which this du- 
ality extends to the reflection factors also. Con- 
sider the triple of coupling constants (/?, do, ax). 



then it is simple to check that the reflection factor 



(3.5) with E and F given by (6.16) and (3.17) 
respectively, is invariant. 

8. Discussion and conclusions 



The sinh-Gordon model is just about the sim- 
plest massive model one might contemplate be- 
yond free field theory. Yet it has been a long 
and tortuous road to discover a dynamical ar- 
gument which would enable us to derive a re- 
lationship between the parameters occurring in 
the Lagrangian formulation of the model and the 
parameters arising from other, purely algebraic 
considerations — in the sense of deductions made 
using the boundary version of the Yang-Baxter 
equations together with the bootstrap. Our ar- 
gument is certainly non-perturbative but it might 
not be exact. We hope that it will turn out to be 
exact but we cannot be sure. 2 

The properties of the boundary breathers we 
have found are quite intriguing. We should like to 
see an explanation for the surprising facts, for ex- 
ample concerning the dependence of the classical 
action on the boundary parameters. We would 
also be interested in exploring the situation with 
two boundaries. There, the field theory is con- 
fined to an interval and the spectrum of states 
within the interval should depend on the bulk 
coupling and four boundary parameters. One 
might expect the energy spectrum to be deter- 
mined by a pair of reflection factors, one for each 
boundary, but that idea requires the assumption 
of factorization. In other words, if the bound- 
aries behave independently we would expect a 
relationship of the type 

e UkL R^(k)R^ L \-k) = 1, (8.1) 

where the two boundaries are situated at x = 
±L, and the corresponding reflection factors are 
R( ±L \ It would be interesting to see how this 



2 It should be remarked that Al. Zamolodchikov has 
also calculated (but not yet published) the parameter de- 
pendence derived using quite different arguments 
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spectrum might be compatible with the quanti- 
zation of periodic classical solutions. 

There are several directions to go beyond 
sinh-Gordon. For example, the next simplest 
model with a single scalar field is the model based 

(2) 

on a 2 data. There, instead of having two bound- 
ary parameters, there is only one but it can arise 
in two distinct ways For this reason, we 

expect a greater variety of boundary breathers. 
The classical reflection factors are known [[ll] but 
at the moment the analogue of Ghoshal's for- 
mula has not been found. Indeed, analysing the 
boundary breathers might be of great assistance 
for this case, to serve as a guide in finding the 
correct expressions. Besides, there is a lot to do 
to classify completely all the reflection data for 
all the affine Toda models, extending what is al- 
ready known for the a! series [[bi] . 
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